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ABSTRACT 

We make a prediction of the cosmic supernova rate history as a composite of the supernova rates in 
spiral and elliptical galaxies. We include the metallicity effect on the evolution of Type la supernova (SN 
la) progenitors, and construct detailed models for the evolutions of spiral and elliptical galaxies in clusters 
and field to meet the latest observational constraints. In the cluster environment, the synthesized cosmic 
star formation rate (SFR) has an excess at z ^ 3 corresponding to the early star burst in ellipticals and 
a shallower slope from the present to the peak at the redshift of z ~ 1.4 compared with Madau's plot. In 
the field environment, we assume that ellipticals form over a wide range of redshifts as 1 ^ z 4. The 
synthesized cosmic SFR has a broad peak around z ~ 3, which is in good agreement with the observed 
one. The resultant cosmic SFRs lead to the following predictions for the cosmic SN la rate: 1) The 
SN la rate in spirals has a break at z ~ 2 due to the low-metallicity inhibition of SNe la, regardless of 
clusters or field. 2) At high redshifts, the SN la rate has a strong peak around z ~ 3 in clusters, whereas 
in field much lower rate is expected, reflecting the difference in the formation epochs of ellipticals. 

Subject headings: cosmology: theory — galaxies: abundances — galaxies: evolution — supernovae: 
general 



1. INTRODUCTION 



The search for high-redshift supernovae has been exten- 
sively conducted mainly to determine cosmo logical param- 
eters by the Supernova Cosmology Project ( Perlmutter ct 
al. 1999|) an d the High-z Supernova Search Team ( Schmidt 
et al. 199§| ). High-redshift supernovae can also provide 



useful information on the star formation history in the 
universe. Type la su pernovae (SNe la) have been discov - 
ered up to z ~ 1.32 ( [GiUiland, Nugent fc PhiUips 19991) 
and the SN la rate has been estimated up to z ~ 0.5 ( |PaiE| 
et al. 19961; |Pain 19991). With the Next Generation Space 



Telescope, both SNe la and Type II supernovae (SNe II) 
will be observed through z ~ 4. The SN II rate is di- 
rectly connected with the star formation rate (SFR), and 
the SN la rate can also trace the SFR combined with a 
SN la progenitor model. In a theoretical approach, the 
cosmic SN la rate as a function of redshift has been con- 
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to z ~ 5 from UV and Ha luminosity densities with the 
help of spectral population synthe sis models (e.g, Madarj 
et al. 1996 ; ponnoUy et al. 1997 ). The observed cosmic 
SI'R by Madau et al. (1996) shows a peak at z ~ 1.4 and 
a sharp decrease to the present. However, UV luminosi- 
ties which are conve rted to the SFRs m ay be affected by 
the dust extinction ( Pettini ct al. 199^ ). Recen t updates 
of the cosmic SFR show some different features (Tresse fc 
Maddox 199^; IGronwall 1998|; |Treyer et al. 1998|; |Hughe; 



et al. 1998| ; |Glazebrook et al. 1999| ; |Steidel et al. 1999| ), 



one of which suggests that a peak lies around z ~ 3. 

Among the several attempts to calculate the cosmic SN 
la rate using the observed cosmic SFR, K98 predicts that 
the cosmic SN la rate drops at z ~ 1 — 2, because the SN 
la occurrence depends on the metallicity. In their SN la 
progenitor model, the accreting white dwarf (WD) blows 
strong winds to reach the Chand rasekhar (Ch) mass limit 
pachisu, Kato fc Nomoto 1996| , 1999b). If the iron abun- 
dance of the progenitors is as low as [Fe/H] ^ — 1, then 
the wind is too weak for SNe la to occur. This model 
successfully reproduces the observed chemical evolution of 
the solar neighborhood such as the evolution of the oxygen 
to iron ratio and the abundance distribution function of 
disk stars (K98). 

Their finding that the occurrence of SNe la depends on 
the metallicity of the progenitor systems implies that the 
SN la rate strongly depends on the history of the star for- 
mation and metal-enrichment. The universe is composed 
of different morphological types of galaxies and therefore 
the cosmic SFR is a sum of the SFRs for different types 
of galaxies. As each morphological type has a unique star 
formation history, we should decompose the cosmic SFR 
into the SFR belonging to each type of galaxy and calcu- 
late the SN la rate for each type of galaxy. 

In this paper, we first construct the detailed evolution 
models for different types of galaxies which are compati- 
ble with the stringent observational constraints, and apply 
them to reproduce the cosmic SFR for two different envi- 
ronments, i.e., clusters and field. Secondly we confirm that 
the metallicity-dependent SN la progenitor model, which 
has been already tested for the solar neighborhood, can 
explain the present supernova rates for all types of galax- 
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ies. Finally combining the above SN la model with the 
self-consistent galaxy models, we calculate the SN la rate 
history for each type of galaxy and predict the cosmic SN 
la rate as a function of redshift. 

In the next section, we describe our computational 
method of the galaxy evolution with our SN la progen- 
itor model. In section 3, we construct the star formation 
histories of spirals and ellipticals, and predict their super- 
nova rate histories in section 4. In section 5, we make a 
prediction of the cosmic supernova rates as a composite of 
different types of galaxies in clusters and field. Discussion 
and conclusions are given in sections 6 and 7, respectively. 

2. MODELS 

2.1. Type la Supernova Model 

The progenitors of the majority of SNe la are most likely 
the Ch mass WDs (e.g., Nomoto, Iwamoto fc Kishimotc] 
1997a; Nomoto et al. 2000| for recent reviews), although 



the sub-Ch mass models might correspond to some pecu- 
liar subluminous SNe la. The early time spectra of the 
majority of SNe la are in excellent agreement with the 
synthetic spectra of the Ch mass models, while the spec- 
tra of the sub-Ch mass models are too blue to be compat- 
ible with ob servations (Hoflich & Khokhlov 1996; Nugent] 
et al. 1997 ). For the evolution of accreting WDs toward 



the Ch mass, two scenarios have been proposed: One is 
a double-degenerate (DD) scenario, i.e., merging of dou- 
ble C-l-0 WDs with a combined mass surpassin g the Ch 



mass limit (Iben fc Tutukov 1984; Webbink 1984), and the 



other is a single-degenerate (SD) scenario, i.e., accretion of 



for a review), but theoretical modeling has indicat ed that 



the merging o f WDs does not make typical SNe la ( Saio 
Nomoto 198^ , 1998; [Segretain, Chabrier fc Mochkovitch| 



1997D , and the lifetime of SNe la progenitors predicted by 



the DD scenario is too short to be compatible with the 
chemical evolution of the solar neighborhood (K98). 

Our SD scenario has two progenitor systems: One 
is a red-giant (RG) companion with the initial mass of 
Mbg,() ^ 1-^0 and an orbital period of tens to hundreds 
days (Hachisu et al. 1996, 1999ab). The other is a near 
main-sequence (MS) companion with an initial mass of 
-^MS. o ^ 2 — 3Aff7i and a period of several tenths of a day 



gions are much smaller for lower metallicity because the 
wind becomes weaker. The wind velocity depends also on 
the luminosity L of the WD. The more massive WD has 
a higher L, thus blowing higher velocity winds (Hachisu 
et al. 1999b). In order for the wind velocity to exceed 
the escape velocity of the WD near the photosphere, the 
WD mass should be larger than a certain critical mass 
for a given [Fe/H]. This implies that the initial mass of 
the WD AfwD,o should already exceed that critical mass 
in order for the WD mass to grow to the Ch mass. This 
critical mass is larger for smaller [Fe/H], reaching I.IMq 
for [Fe/H] = -1.1 (Fig.l of K98). Here we should note 
that the relative number of WDs with A fwn.n ^ l.lMp, i s 



quite small in close binary systems (Umcda et al. 1999). 

> 1. 2M(7i, the accretion lead s to collapse 

Therefore, 
-1.1 in our 



WD,0 



And for Afi 

rather than SNe la ( [Nomoto fc Kondo 1991D 
we assume that no SN la occur at [Fc/H] < 
chemical evolution model. 

The lifetimes of SNe la (i.e., the lifetimes of the binary 
systems from the formation through the explosion) are 
determined from the main-sequence lifetimes of the com- 
panion stars with the initial masses in the range rrid^e < 
m < md,ti. In our chemical evolution model, we adopt 
-^^WD,o — lAf© and Z = 0.004 as a representative case. 
ThenTOd,£ = 0.9AfQ and TOd,« = l.SAf© for the RG-I-WD 
system and m^j = 1.8Af0 and md.u = 2.6AfQ for the 
MS+WD system, as seen in Figure 1 of K98. The distri- 
bution function of the companion stars (mass donors) is 
assumed to be 0d(w) oc m^^ with the slope of a; = 0.35, 
and normalized to unity at Wd.^ < 'tt, < md.w This slope 
is determined from the distribution function of the initial 
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mass ratio distribution in binaries (Duguennoy fc Mayor 



The progenitors of SNe la are assumed to have main 
sequence masses in the range of rrip^i = 3Af0 and nip^u = 
8Mq and form C-l-0 WDs. The total fraction of primary 
stars which eventually produce SNe la is denoted by b 
(see Eq.|l^ in section 2.2 for the exact definition of b). 
The parameter b may depend on the binary fractions. If 
the binary fraction depends on the gas density in the star 
forming region, there should be a difference in b between 
spirals and ellipticals (we can expect the higher fraction 
for ellipticals). However the number of low-mass X-ray 
binaries per unit blue luminosity appea rs to be almost 



to se^iral days (pji fc van den Heuvel 1997 ; Hachisu et 
ah 1999a In tnese progenitor systems, the C-l-O WD 



the same between spiral s and ellipticals (Canizares, Fab- 
biano fc Trinchieri 1987|; [Matsushita et al. 1994] ), which 



accretes H-rich matter from the companion star at a suf- 
ficiently high rate (10^^ — 10~^A^o yr~^) to increase its 
mass to the Ch mass through steady hydrogen burning. 
With such a rapid accretion, the WD blows optically thick 
strong winds that play a key role in stabilizing the mass 
transfer and avoiding a common envelope formation. 

The optically thick winds are driven b y a strong peak of 
OPA L opacity at logT(K) ~ 5.2 (e.g., Iglesias fc Rogers 
1993[ ). Since the opacity peak is due to iron lines, the 
wind velocity v^, depends on the iron abundance [Fe/H] 
(K98; [Hachisu fc Kato 200(]| ), i.e., is higher for larger 
[Fe/H]. The metallicity effect on SNe la is clearly demon- 
strated by the size of the regions to produce SNe la in the 
diagram of the initial orbital period versus initial mass of 
the companion star (see Fig. 2 of K98). The SN la re- 



means that the binary fraction does not depend on the 
gas density in the star-forming region. The parameter b 
may also depend on the metallicity. We have not included 
the metallicity dependence of the size of the SN la region. 
According to the simulations of Hachisu fc Kato (2000), 
the area in the orbital period-companion mass plane does 
not increase with the metallicity for the MS-I-WD systems, 
but does increase for the RG+WD systems. Thus, only b 
for the RG-I-WD systems 6rg may be larger for higher 
metallicity. However there is no significant difference be- 
tween the mean iron abundances of progenitors for the 
RG-I-WD systems in ellipticals and spirals. Thus, we de- 
termine the values of b to reproduce the chemical evolution 
in the solar neighborhood, and adopt the best-fit values for 
all types of galaxies. (We also show the results for other b 
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parameters in Appendix.) We search the parameter b with 
the test for the [0/Fe]-[Fe/H] relation and the abun- 
dance distribution in the solar neighborhood, and get the 
allowed values of 0.04 < &ms - 0.055 for the MS+WD sys- 
tems and 0.01 < &rg £ 0.04 for the RG+WD systems un- 
der 2&MS + ^RG ~ 0.12. The best-fit values are 5ms = 0-05 
and 6r,g — 0.02 with the probability of more than ~ 80%; 
the probabihty is as low as ~ 65% for 6ms = 5rg = 0.04^ 
(see Appendix). The fraction of primary stars which even- 
tually produce SNe la is given as functions of the compan- 
ion's mass md and the companion's lifetime (i.e., lifetime 
of SNe la progenitor tia) as 



1 r"'^.^ 1 

Bm{md) = b 0d(m-d) W — 4>d{m) dm 



rrid 



(1) 
(2) 



Figure |l| shows Bm{md) (upper panel) and Bt{tia) (lower — 
panel) for the RG-hWD and the MS-hWD systems. In the 
lower panel, the solid and the dashed lines show the cases 
for Z = 0.002 and Z = 0.02 (solar), respectively. 

In the one-zone uniform model for the chemical evolu- 
tion of the solar neighborhood, the heavy elements in the 
metal-poor stars originate from the mixture of the SN II 
ejecta of various progenitor masses. The abundances aver- 
aged over th e progenitor masses of SNe II predicts [0/Fe ] 
~ 0.45 (e.g., [Tsujimoto ct al. 1995| ; [Nomoto ct al. 1997bD . 
Later SNe la start ejecting mostly Fe, so that [0/FeJ dc- 
creases to ~ around [Fe/H] ~ 0. The low-metallicity 
inhibition of SNe la predicts that the decrease in [0/Fe] 
starts at [Fe/H] '--^ — 1. Such an evolution of [0/Fe] well 
explains the observations (K98). 

We should note that some anomalous stars have [0/Fe] 
'--^ at [Fe/H] — 1. The presence of such stars, however, 
is not in conflict with our SNe la models, but can be under- 
stood as follows: The formation of such anomalous stars 
(and the diversity of [0/Fe] in general) indicates that the 
interstellar materials were not uniformly mixed but con- 
taminated by only a few SNe II (or even single SN II) 
ejecta. This is because the timescale of mixing was longer 
than the time difference between the supernova event and 
the next generation star formation. The iron and oxygen 
abundances produced by a single SN II vary depending on 
the mass, energy, mass cut, and metallicity of the progeni- 
tor. Relatively smaller mass SNe II (13— I5M0) and higher 
explo sion e nergies tend to produ ce [O /Fe] ~ ( Nomoto et 
al. lj97b|; lUmeda et al. 20"00| ). Those metal poor stars 
with [0/Fe] ^ may be born from the interstellar medium 
polluted by such SNe II. 

The metallicity effect on SNe la can also be checked 
with the metallicity of the host galaxies of nearby SNe la. 
There has been no evidence that SNe la have occurred 
in galaxies with a metallicity of [Fe/H] ^ — 1, although 
host galaxies are detected only for one third of SNe la and 
the estimated metallicities of host galaxies are uncertain. 
Three SNe la are observed in low-metallicity dwarf galax- 
ies; SN1895B and SN1972E in NGC 5253, and SN1937C in 
IC 4182. Metallicities of these galaxies ar e estimated to b e 
[0/H] = -0.25 and -0.35, respectively ( [Kochanek 1997| ). 

^In K98, we adopt 6ms = ''RG = 0.04 which is determined from the 



If [0/Fe] - as in the Magellanic Clouds, [Fe/H] 0.25 

and —0.35 which are not so small. Even if these galaxies 
have extremely SN II like abundance as [0/Fe] ^ 0.45, 
[Fc/H] ^ —0.7 and —0.8 (being higher than —1), respec- 
tively. Since these host galaxies are blue {B — V = 0.44 for 
NGC 5253 and S-y = 0.37 for IC 4182 according to RC3 
catalog), the MS+WD systems are dominant progenitors 
for the present SNe la. The rate of SNe la originated from 
the MS-I-WD systems is not so sensitive to the m etallicity 
as far as [Fe/H] > -1 ( [Hachisu fc Kato 20001 ). Even if 
[Fe/H] ~ —0.7 in such blue galaxies, therefore, the SN la 
rate is predicted to be similar to those in more metal-rich 
galaxies. 

Futhermore, our SN la model can well explain the di- 
versity of SN la luminosity, by assuming that the more 
luminous SNe la are produced from smaller M-wd,o be- 
cause of the larger C/0 ratio. Observationally, the most 
luminous SNe la occur only in spirals, whil e both spirals 



and e llipticals are hosts for dimmer SNe la (Umeda et al 



1999). In our progenitor model, the WD with small A/wd,o 



(~ O.7M0) can reach the Ch mass only when the compan- 
ion stars is massive enough (i.e., iia is short enough) to 
supply the necessary amount of mass (~ O.7M0). Thus 
the brightest SNe la can occur only near the star forming 
region. The variation of the peak brightness of SNe la has 
been found to depend also on the location in the galaxy, 
i.e., the diversity decreases toward the ou ter regions of the 
galaxy (Wang, Hoflich & Wheeler 1997). In our model, 
this trend is due to the metallicity gradient of the galaxy. 
In the outer region, the metallicity is so small that Mwd.o 
should be large enough to blow strong winds; thus only 
dim SNe la appear. 

2.2. Galactic Evolution Model 

A simplified model of the galactic c hemical evolu tion is 
constructed with analytical equations ( Tinsley 1980 ). The 
accretion of primordial gas from a reservoir in a halo is as- 
sumed to form a galaxy. By assuming that the accretion 
per time t (the infall rate i?in) is proportional to the mass 
of the reservoir, we get the exponential form of the infall 
rate with a infall timescale Ti as 



i?i„ = — exp(-t/Ti). 

Ti 



(3) 



Here we define the total baryon mass of the galaxy as the 
mass of the reservoir. 

The gas fraction /g is defined as the ratio of the gas 
mass in the bulge and disk to the total baryon mass of the 
galaxy (bulge, disk and halo), and Zi is the mass fraction 
of a heavy element i in the gas. For ellipticals, the gas 
includes the gas ejected in the galactic wind. The_SFR__0 
is ass umed to be proportional to the gas fraction ( ^chmidt 
T959I ) as 

where Tg is the star formation timescale. Then the time 
variations of /g and Zi are given by the following equa- 
tions: 



dt 



-ij + E + Eu 



Ri, 



(5) 



Kolmogorov-Smirnov (KS) test. 
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dt 



I g,, I E.^ji I E,a, I Z,,m-Rm.(6) 



range fro m nii — 0.05Mq to = 50Mq (Tsujimoto ei 



The metallicity Z is defined as the sum of from C to 
Zn. The initial conditions are /g = and Zi = 0. The 
metalhcity of the infaU gas Zi^in is assumed to be 0. 

From dying stars, gas is ejected into the interstellar 
medium by mass loss and SNe II at a rate of E and by 
SNe la at a rate of Ei^. Heavy elements are ejected at a 
rate of E^^, Ez^^n, and Ez^^ia by mass loss, SNe II, and 
SNe la, respectively. These ejection rates are given by the 
following equations; 

E^ {I - Wm) i^it - Tm) 4>{rn) dm, (7) 



Ez^^ I {l-Wrn-Pzirn,ll) Z^{t-Trn)i^{t-Trn)4>{m) dm, 
J nit 

(8) 

Ezi,u^ / Pz,m,uij{t - Tm)(l){m) dm, (9) 

J nit 

-Bia = mco 7?.Ia, (10) 

£^z.,Ia = "^COP^im,Ia7?.Ia■ (H) 

Here mco = 1.38Mq is the white dwarf mass at the ex- 
plosion. TZia, denotes the SN la rate, which is obtained 
as 

T^ia — b — (?!)(m) dm 

"'inaxfmp^^, mt] ^ 

— ipit ~ Tra)4' dim) dm. (12) 

max[md.^ , mt] ^ 

As noted in section 2.1, our SN la scenario has two types 
of progenitors (i.e., the MS+WD and the RG+WD sys- 
tems). We calculate the SN la rate for each binary with 
each b, m^^e, and m^.u, and combine them. The SN II rate 
TZii is also obtained as 

p. 1 
TZii= / — '4>{t ~ Tm) 4>{m) dm. (13) 

The lower mass limit for integrals is the turning off mass 
TOj at t which is the mass of the star with the main se- 
quence lifetime Tm — t. Tm is taken from Kodama & 
Arimoto (1997) as a function of metallicity Z. Wm is the 
remnant mass fraction, which is the mass fraction of a neu- 
tron star or a white dwarf. Pzim.ii and Pzim.ia are the stel- 
lar yields which are the mass fractions of newly produced 
and ejected heavy element i, whic h are given from the su 
pernovae nucleosynth esis model ( [Tsujimoto et al. 1995 
Nomoto ct al. 1997bD with Pz.m.ii = for m < lOM©. We 



do not include the dependence of u;,„, Pzim,ii and Pz-m.ia 
on the stellar metallicity. 

The initial mass function (IMF) is assumed to have time- 
invariant mass spectrum (j){m) cx normalized to unity 
at mi < m < mu. Theoretical arguments indicate that 
the IMF originates from fragmentation of a gas cl oud al- 
most independent l y of local physics in the gas ( [Low & 
Lynd 4n-BeU 1976| ; |Silk 1977| ). A solar-neighborhood IMF 
would therefore be a good ap proximation, a nd we adopt 
the Salpeter slope of x = 1.35 ( ^alpeter 1955 ) and a mass 



al. 1997) 



The time variation of the stellar fraction /g is calculated 



as 



dt 



= i,-E-Eia, 



(14) 



and the mean stellar metallicity Z^.s at the time t is ob- 
tained from the conservation of heavy elements; 

Z^ /g + Zt., fs^ f {Ez^Jl + Ez,M + Z^,in Rin} dt. (15) 

Jo 

The photometric evolution of galaxies is calculated from 
the summation of the simple stellar population, which is 
defined as a single generation of coeval and chemically ho- 
mogeneous stars of various masses, and taken from Ko- 
dama & Arimoto (1997) as a function of age t and metal- 
licity Z. The passbands of photometric systems and the 
zero points are the same as Kodama & Arimoto (1997). 

For a standard model, we adopt Hq = 50 km s^^ 
Mpc~^, r^o = 0.2, Ao = 0, and the galactic age of 15 
Gyr, which corresponds to the redshift at the formation 
epoch of galaxies of Z{ ^ 4.5. 

3. STAR FORMATION HISTORIES OF GALAXIES 

3.1. Spiral Galaxies 

The star formation history can be inferred from the ob- 
served present-day colors of galaxies, by using the well- 



known technique of stell ar population synthesis ( Arimoto 
Yoshii fc Takahara 1992 ). In this paper, we determine the 
timescales Tj and Tg in equations (^ and (^) to reproduce 
both observed colors and gas fractions for four types of 
spirals, i.e., SOa-Sa, Sab-Sb, Sbc-Sc, and Scd-Sd. These 
values are summarized in Table |^. As shown in the top 
panel of Figure H, earlier types of spirals form larger frac- 
tions of stars at an early epoch, thereby having redder and 
smaller gas fractions at present. The excellent agreements 
between models and the observations are shown in mid- 
dle and bottom panels of Figure ^. The resultant present 
metallicities and colors are tabulated in Table ||. Observa- 
tional data of the present B — V colors for various types 
of spirals are taken from Roberts & Haynes (1994). We 
use the gas (i.e., HI-I-H2) fractions which are normalized 
by the present blue luminosity of the galaxy to avoid the 
uncertainty in the fractions of the dark matter. The HI 
mass is taken from Roberts & Haynes (1994), and the H ^ 
mass is derived from the H2/HI ratios ( Casoli et al. 1998 ). 



3.2. Elliptical Galaxies 

The star formation history in elliptical galaxies is still 
controversial, with single star burst models on one hand 
and continuous star formation models at the other ex- 
treme. In the former, elliptical galaxies are formed through 
dissipative collapse of a protogalactic c loud with a s i ngle 



star burst at a very early epoch (e.g., Larson 1974 ; Ari 



moto k Yoshii 1987| ; [Kodama fc Arimoto 19971 ). In the 



latter, elliptical galaxies grow through mergers of gaseous 
galaxies with a continuous star formation through later 
epoch (e.g. 
1998[). 



Kauffmann & Chariot 1998; Baugh et al 



The dissipative collapse scenario assumes that the star 
formation in ellipticals has stopped by a loss of gases due to 
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a sup ernova -driven galactic wind (Larson 1974; Arimoto & 
Yoshi 



1987), so that the bulk of the stars are old and have 
formed at z 2. The galactic wind model can well repro- 
duce the passive evolution of colors observed in cluster el- 
liptic als (Stanford, Eisenhardt fc Dickinson 1998| ; Kodama| 
& Ari moto 1997| ). On the other hand , the semi-analytical 



simulati ons of hierarchical c lustering ( Kauffmann fc Char _ 
lot 19p8; Baugh et al. 1998| ) based on the cold dark matter 
(CDM) scenario suggest that the formation of ellipticals 
are protracted toward lower redshifts especially in the low- 
density region such as present-day field. Such environmen- 
tal effect is imprinted in the morph ology-density relation 
(Dressier 1980; Dressier et al. 1997). There are some ob- 
servational evidences that not al l field ellipticals are as old 



as cluster ellipticals (Zepf 1997 ; Franceschini et al. 1998 
Kodama. Bower k, Bell 1999] ). However the relation of 



Mg2 indices and velocity dispersions implies the age dif- 
ference of only ^ 1 Gyr amon g cluster, group and field 
ellipticals (Bernardi et al. 1998). 

Taking into account the possibility that there exists the 
environmental effect on the galaxy formation, we set two 
models for ellipticals in clusters and field, respectively. 
1) For cluster ellipticals, we adopt the galactic wind model. 
The epoch of galactic wind tg^ (i.e., the epoch of the end 
of star formation) is determined from the dynamical po- 
tential of the galaxy (Larson 1974), and we assume igw — 1 
Gyr on the average which corresponds to the redshift of 
z ^ 3. The adopted Tj and Ts in equations (||) and (^) are 
summarized in Table |l|, and the resultant present metallic- 
ities and colors are in Table This star formation history 
is constructed to reproduce the observational constraints 
such as the present mean stellar m etallicity of ^ O.GZp 



averaged over the whole galax y (Kobayashi fc Arimotj 
1999|J the present B ~V color ( [Roberts fc Haynes 199^ , 



and the color evolution observed at ^ z 1 as shown in 
Figure ||a. The circles, triangles, and crosses show the av- 
erage differences in the observed colors of cluster galaxies 
at each redshifts relati ve to the same rest-frame colors of 
Coma E-t-SO galaxies ( |Stanford, Eisenhardt fc Dickinso"ij 
19980 Thus, no color difference represents no evolution 
relative to Coma. The observed color differences are as 
small as our model prediction, which means that ellipti- 
cals have evolved passively from z ^ 1. In this model, the 
duration of the star formation is short enough to be con- 
sistent with the observational estimates of the magnesium 
to iron ratio [Mg/Fe] ~ 0.3 (Worthey, Faver fc Gonzalez 



1992 1; [Kobayashi fc Arimoto 1999D " 
2) For field ellipticals, we adopt the same star formation 
history as in cluster ellipticals, and assume that the for- 
mation epochs span a wide range of redshifts. The dis- 
tribution function of the formation epoch (j)zf is assumed 
as 



2^/3 
V3- 1 



2a/3 



(16) 

which gives (f>zf = at z = and = at z = 2, and 
is normalized to unity at < z < Zf. As shown in Figure 
1^, this formulation can meet the distribution function of 
Zf which is derived from the observational estimates of Zf 
for 34 field elliptical s in the Hubble Deep Fie ld, using the 
broadband spectra (Franceschini et al. 1998). 



4. SUPERNOVA RATES IN GALAXIES 

Present s upernova rates observed in the various types 
of galaxies ( [Cappellaro et al. 1997 ; Cappellaro, Evans & 
Turatto 1999) put the constraints on the SN la progenitor 
models. Using the galaxy models in section 3, we show 
that our SN la model can well reproduce the present su- 
pernova rates in both spirals and ellipticals. 

4.1. Spiral Galaxies 

The observed SN H rate TZu in late-type spirals is about 
twice the rate in early-type spirals. On the other hand, 
the observed SN la rates TZu in both types of spirals are 
nearly the same. Thus the present Tiia_/Tiu ratio in early- 
type spirals is larger than that in late-type spirals by a 
factor of '-^ 2. Such a difference in the relative frequency is 
a result of the difference in the SFR (see section 3.1), be- 
cause the different dependences of TZu and TZin on the SFR 
are due to the different lifetimes of supernova progenitors. 
Therefore the observed TZii^/TZu ratio gives a constraint on 
the SN la progenitor model. 

Figure shows the predicted evolutionary change in 
Ti-ia/TZii ratios for early and late types of spirals, com- 
pared with the observations. The solid, dashed, and dot- 
ted lines show the results for our SN la model, the single 
delay-time model with tig_'^ 1.5 Gy r (Yoshii, Tsujimoto fc 
Nom oto 1996), and the DD model (Tutukov fc Yimgelson 
1994 ), respectively. In the DD model, the lifetime of ma- 
jority of SNe la is ^ 0.1 — 0.3 Gyr, so that the evolution of 
TZia, is similar to TZu. Therefore TZis,/TZii is insensitive to 
the SFR. This results in the small differences in TZin/TZu 
among the various type of spirals, which is not consistent 
with observations. For the similar reason, the single delay- 
time model with tja ~ 1.5 Gyr is not acceptable. 

In our SN la model, if the iron abundance of progenitors 
is [Fe/H] }t — 1, the occurrence of SNe la is determined 
from the lifetime of the companions, which is iia ~ 0.5—1.5 
Gyr for the MS companions and ^2 — 20 Gyr for the RG 
companions (see Fig.|l|). If SNe la occurred only in the 
MS+WD systems with relatively short lifetimes, TZin/TZu 
would have been insensitive to the SFR. On the contrary, if 
SNe la occurred only in the RG-f WD systems with longer 
lifetimes, the present difference in TZia/TZu between early 
and late type spirals would have been too large, reflecting 
the large difference in SFR at an early epoch. Owing to 
the presence of these two kinds of the progenitor systems 
in our SN la progenitor model, the observed difference in 
TZia/TZu can be reproduced. We note that the present SN 
la rate depends on the b parameters, especially &rg (see 
Appendix for more details). 

4.2. Elliptical Galaxies 

Figure ^ shows the SN la rate history in ellipticals. As 
noted in section 3.2, we assume that the bulk of stars 
in cluster ellipticals are formed at z 3 and have ages 
older than 10 Gyr. Thus, in the single delay-time model 
(dashed line) and the DD model (dotted line), the SN la 
lifetimes are too short to reproduce the observed SN la 
rate at the present epoch. In ellipticals, the chemical en- 
richment takes place so early (see Fig.|^) that the metal- 
licity effect on SN la is not effective, and therefore the 
SN la rate depends almost only on the lifetime. Our SN 
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la model (solid line) includes the RG+WD systems with 
ha ^ 10 Gyr, thus well reproducing the present SN la rate 
in ellipticals. 

A burst of SNe la occur after ^ 0.5 Gyr from the begin- 



(19) 



ning elf the HLai' foi'maLioii, becauae SNe la atait Lo occur 
from the MS+WD sysLems aL Llieii shuilesL lifetime. The 
second peak of the SN la rate appears after ~ 2 Gyr, be- 
cause of the onset of SNe la from the RG+WD systems. 
If we apply the age-redshift relation in Figure ^, these two 
peaks appear at z 2.5 and ~ 1.5, respectively. From 
z ~ 0.2 to z '--^ 0, the SN la rate gradually decreases by 
- 40%. Majority of SN la progenitors with Z = 0.002 
have already exploded at z < 0.2 and only more metal- 
rich SNe la occur at z ~ 0. This is because the lifetime 
r,„ of the smallest mass companions (0.9Mq) depends on 
the metallicity as ~ 11 Gyr for Z ~ 0.002 and ~ 19 
Gyr for Z = 0.02 (see Fig.Q). 

The decrease in the SN la rate from z ^ 0.2 to z ^ 
depends also on the star formation history in ellipticals 
and the galactic age. If we adopt the galactic age of 12 
Gyr, such a decrease in the SN la rate does not appear. 
If ellipticals have undergone the relatively continuous star 
formation, as suggested by the hierarchical clustering sim- 
ulations, the SN la rate might keep constant to the present. 
The predicted rate for Hq = 50 km s^^ Mpc^^ adopted in 
our calculation is a little higher than the observations at 
present. We can get better fit for Hq — 65 km s^^ Mpc^^, 
because the absolute value in our model is independent 
of the cosmological parameters, while the observed SN la 
rate is proportional to Hq. The present SN la rate also 
depends on the b parameter, especially 5rg (see Appendix 
for more details). 

5. COSMIC SUPERNOVA RATE 

Galaxies that are responsible for the cosmic SFR have 
different timescales for the heavy-element enrichment, and 
the occurrence of supernovae depends on the metallicity 
therein. Therefore we calculate the cosmic supernova rate 
by summing up the supernova rates in spirals and ellipti- 
cals with the ratio of the relative mass contribution. The 
relative mass contribution of z-th type of galaxies is ob- 
tained from the observed relative lumin osity propor tion 
d = 0.215,0.185,0.160,0.275, and 0.165 ( pence 1976| ) for 
ellipticals, SOa-Sa, Sab-Sb, Sbc-Sc, and Scd-Sd, respec- 
tively, and the calculated mass to light ratio in B-band 
{M/LB)i for each galaxy model, as given in Table I 

To compare the observed cosmic SFR, we convert the 
predicted cosmic SFR per mass to the cosmic SFR per 
volume by multiplying a constant (pc^^gcx)) as follows; 



V'c 



[Moyr-iMpc-''] = p,f7gc 



S, c, [M/Lsh [Gyr- 

T.,C,{M|LB)^ 



(17) 

Here pc is the critical density in the cosmology (pc = 
^qHq) and Ogoo is the initial comoving density of the 
gas defined by Pei & Fall (1995). figoo is constrained from 
the stellar comoving density at present. The gas and stel- 
lar comoving densities are given by 



Si Ci {M/LB)i fgA 



We adopt Vl^oo — 3.5 x 10~^ft.^^ to repr oduce the present 



stellar comoving density logSlso — —2.3 (Fukugita, Hogan 



& Peebles 1998), which gives the present gas comoving 
density of logfigo = —3.6. 

The cosmic SN II or la rate T^sn per luminosity and per 
volume are given as 



.cosmic 

[(lOi°Ls0)"^Century" 



T^SN, cosmic [yr ^MpC ^] = Pc^^g 



-ij ^ S,c,7^sN,^ [(10i"£BQ)-iCei 

Si Cj 

(20) 

c, {M/Lb\ 7^sN,^ [Mq iGyr-1] 



Si a {M/LB)i 
(21) 

If we adopt larger figoo , we get larger fiso and thus larger 
V'cosmic- However, the 7?.sn. cosmic per luminosity does not 
depend on figoo. 

5.1. In Clusters 

First, we make a prediction of the cosmic supernova 
rates in cluster galaxies. We use the galaxy models which 
are in good agreements with the observational constraints 
as shown in Figures I and |. Figure |^a shows the cos- 
mic SFR along redshift (solid line) as a composite of those 
in spirals (long-dashed line) and ellipticals (short-dashed 
line). In our galaxy models, ellipticals undergo a star burst 
at z 3 and the duration of the star formation is ^ 1 Gyr, 
while spirals undergo relatively continuous star formation. 
Thus, only the SFR in spirals is responsible for the cosmic 
SFR at z ;$ 2. For comparison, the observed cosmic SFR 
in field, so-called Madau's plot ( Gallego et al. 199£ Lilly 
et al. 1996| [Madau et al. 1996| ; ponnoUy et al. 1997] ), is 
also plotted. Compared with Madau's plot, the predicted 
slope from the pre sent to the peak at z ^ 1.4 is a little 
shallower (see also Totani, Yoshii &: Sato 1997 ). Also the 
predicted peak of the SFR in ellipticals at z 3 is not seen 
in Madau's plot. 

Possible sources of these differences between the theo- 
retical SFR and observational data might be as follows; 1) 
Observational uncertainties: Recent observations have re- 
vealed the controversial situation, i.e., the slope of the SFR 
from z ~ to z ~ 1 is much shallower than the Madau's 
plot according to Tresse & Maddox (1998), Treyer et al. 
(1998), and Gronwall (1998), while it is steeper accord- 
ing to Rowan-Robinson et al. (1997) and Glazebrook et 
al. (1999). Also the SFR is higher than Madau's plot at 



2 (Bughes et al. 199S; Pettini et al. 1998) and con 



goo 



(M/L 



B]i 



(18) 



tinues to be high toward z ^ 5 ( |Steidel et al. 1999 ). 2) 
The star fo rmation at z 3 m ay be hidden by the dust 
extinction (Pettini et al. 1998). 3) Ellipticals may have 
formed at z ~ 5 (Totani et al. 1997). We here propose one 
more possibility: 4) The theoretical model is constructed 
for cluster galaxies, while the observational data are taken 
from field galaxies. Then these rates can be different, if 
the star formation histories depends on their environment. 
The cosmic SFR histories in field is much less clear than 
in clusters. 

Figure ^ shows the cosmic supernova rates (solid line) 
as a composite of those in spirals (long-dashed line) and 
ellipticals (short-dashed line). The upper and lower three 
lines show the SN II and la rates, respectively. Our SN 
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la mo del is in Rood agreement with the observed rate at tini ct al. 199^ ). The corresponding cosmic chemical evo 



0.5 ( Pain et al. 1996| ; |Pain 1999| ); here the number 



of sample has been increased from 3 SNe la in Pain et al. 
(1994) to 38 SNe la in Pain (1999), so that the errorbar 
has been greatly reduced. As we will observe supernovae 
at higher and higher redshift, the SN la rates are predicted 
to sharply decrease at certain redshifts, which we denote 
as 2;ia,si ^ia,cj Eind zia,tot for in spirals, ellipticals, and all 
types of galaxies (spirals plus ellipticals), respectively. 

The SN la rate in spirals drops at zia.s ~ 1-9 because 
of the low-metallicity inhibition of SNe la. We can test 
the metallicity effect by finding this drop of the SN la 
rate in spirals, if high-redshift SNe la at z ^ 1.5 and their 
host galaxies are observed with the Next Generation Space 
Telescope. 

In ellipticals, the chemical enrichment takes place so 
early that the metallicity is large enough to produce SNe 
la at z ^ 2. Since we assume that the formation epoch 
of ellipticals is Zf ~ 4.5, the SN la rate in ellipticals de- 
creases at zia,c ^ 2.6. This redshift is determined from 
the shortest lifetime of SNe la, which is ~ 0.5 Gyr in our 
model (Fig||). We should note that zja^o depends on the 
formation epoch of ellipticals, which is constrained from 
the color-mag nitude relation in the clust er elliptical galax- 



ies at z 3 (Kodama & Arimoto 1997). The two peaks 



of SN la rates at z ~ 2.5 and z ~ 1.5 come from the 
MS+WD and the RG-f WD systems, respectively (see sec- 
tion 4.2). Between these peaks, the SN la rate decreases 
at zia,o ~ 1-6. (Note, however, this drop disappears if we 
adopt Zi ^ 3, because the peak from the MS+WD systems 
moves to lower redshifts.) At z 2.5, the total SN la rate 
decreases as ellipticals at zia.tot ^ 2.6. At z ^ 2, the total 
SN la rate decreases as spirals at zia,tot 1-9, because the 
contribution of ellipticals in the total SN la rate is smaller 
than spirals. However, the contribution of ellipticals can 
be larger depending on b (see Appendix). Figure^ is the 
same as Figure ^d, but for the supernova rate per volume. 

These redshifts, zia,s, zii^,c and zia,tot, depend on the 
speed of the chemical enrichment in the host galaxies. 
Figure ^ shows the evolution of the iron abundance of 
gases in spirals (long-dashed line) and ellipticals (short- 
dashed hne). For ellipticals, the model hue is truncated at 
z = 3.05, where ellipticals lose their gases by the galactic 
wind. As the chemical enrichment in ellipticals takes place 
in short time, the iron abundance reaches [Fe/H] '--^ — 1 at 
z ~ 4.3. Then, the low metallicity effect on SNe la does 
not appear. In spirals, the iron abundance reaches [Fe/H] 
~ — 1 later at z ~ 2.3, which results in the drop of the SN 
la rate at zia,s ^ 1-9. 

For comparison, we calculate the cosmic supernova rates 
for the cosmic SFR constructed by connecting the the ob- 
servational data (Figj^a). In Figure the solid line (and 
the dotted line) is obtained with (and withou t) taking into 
account a correction of the dust extinction ( Pettini ct alj 
19980 The initial gas comoving density is assumed to be 
f^goo = 3.2 X lO^^h^^. The adopted cosmic SFR increase s 
from the Ha data point at z = ( Gallego et al. 1995| ), 
makes a peak at z 1.4, and decreases towards higher 
redshifts. It goes through the HDF UV data points which 
take into account a correction of the dust extinction ( Pet 



lution is shown by the dotted line in Figure [7|a. Because of 
the lower SFR than in Figure |^ at high redshifts, [Fe/H] 
~ — 1 is reached at smaller redshift and the total SN la 
rate decreases at zja.tot ^ 1.6^. 

5.2. In the Field 

We also predict the cosmic supernova rates assuming 
that the formation of ellipticals in fields took place for 
over the wide range of redshifts (FigJ|b). Figure ga shows 
the cosmic SFR along redshift (solid line) as a composite 
of those in spirals (long-dashed line) and field ellipticals 
(short-dashed line). The SFR in spirals is the same as Fig- 
ure ^a, but the star formation in ellipticals continues up 
to the present. Although the timescale of star formation 
in each elliptical is as short as in Figure the SFR av- 
eraged over the field ellipticals gradually increases toward 
higher redshifts. This is because the formation epochs of 
the field ellipticals distribute as in Figure |b. The synthe- 
sized cosmic SFR is in good agreement with t he observed 
one, except for the recent Ha data at z ^ 0.9 ( |Glazebrook 
et al. 1999). The peak of the star formation appears at 



3, which is consistent with the recent sub- mm data 



( [Hughes et al. 1998| ). 

Figure ^ shows the cosmic supernova rates (solid line) 
as a composite of those in spirals (long-dashed line) and 
ellipticals (short-dashed line). The upper and lower three 
lines show the SN II and la rates, respectively. As in Fig- 
ure ||b, the SN la rate in spirals drops at zja.s ~ 1-9. The 
averaged SN la rate in ellipticals decreases at zia,o ~ 2.2 
as a result of ~ 0.5 Gyr delay of the decrease in the SFR 
at z ^ 3. (We should note that zia,o depends on the dis- 
tribution function of the formation epochs in Figure ^o, 
which is still uncertain.) Then, the total SN la rate de- 
creases gradually from zia,tot ~ 2 to z ~ 3. However, the 
contribution of ellipticals can be larger depending on b (see 
Appendix). Figure ^ is the same as Figure ||b but for the 
supernova rate per volume. 

The rate of SNe II in ellipticals evolves following the 
SFR without time delay. Then, it is possible to observe 
SNe II in ellipticals around z ~ 1. The difference in the 
SN II and la rates between cluster and field ellipticals re- 
flects the difference in the galaxy formation histories in the 
different environments. 

6. DISCUSSION 

We have made predictions in Figures |^b and |^b that the 
SN la rates have breaks at zia,s, zia,ej and zia,tot in spirals, 
ellipticals, and all types of galaxies, respectively. The ex- 
act values of the break redshifts depend on still uncertain 
factors. Here we discuss these dependences and uncertain- 
ties in order to know what we can learn when the Next 
Generation Space Telescope will obtain these redshifts. 



In the model in Figure 4 of K98, zja 
and Ogoo = 2 X 10~^h~^ was adopted. 



6.1. SN la Rates in Spiral Galaxies 

For spirals, zja.s corresponds to the point where the 
metallicity reaches [Fe/H] ~ — 1. Therefore zia,s is deter- 
mined by the Z — z relation (metallicity-redshift relation), 
which involves the uncertainties in the t—Z relation (age- 
metallicity relation) and the t — z relation (time-redshift 

is as small as ~ 1.2, because the dust correction was not included (like the dotted line in Fig.pk) 
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relation). Note the t — z relation in galaxies depends not 
only on the cosmological model but also on the formation 
epoch of the galaxies Zf. 

Among the above factors, the uncertainty in modeling 
the evolution of spiral galaxies is relatively small. The 
standard evolution model for spiral galaxies has been ob- 
tained by constructing the detailed model for the Galaxy 
which has a number of observational constraints. Based 
on the parameters on the IMF and the SN la rate deter- 
mined for the Galaxy, we have developed the models for 
other types of spiral galaxies. Futhermore, the parameters 
on the SFR have been constrained from the gas fractions 
and colors of present-day spirals as in Figure 0. 

As for the t — z relation, a set of the cosmological 
parameters (Hq, ^o, Xq) and Zf are constrained from 
the Galactic age which settles dow n to ~ 12 — 15 Gyr 
from the recent HIPP ARCOS data ( pratton et al. 1997 ; 
Harris et al. 1997). (Here we assume the same Zf 
for all types of spirals.) In our calculations, we adopt 
(ifo,ilo,Ao) = (50,0.2,0) that gives the cosmological 
age of 16.5 Gyr and thus allows the galactic age of 
igai = 15 Gyr; these parameters lead to zja.s ~ 1.9. If 
we adopt (i7o,r2o,Ao) = (50,1.0,0) and (65,o'.2,0), which 
give shorter cosmological ages and allows the galactic age 
of iago — 12 Gyr, we obtain smaller values of zia,s ~ 1-3 
and 1.8, respectively. On the other hand, the A-dominant 
model, e.g., (i?o,^^o,Ao) = (65,0.2,0.8) with the galactic 
age of tagc — 15 Gyr, predicts zia,s ~ 2.1. The uncertainty 
in the t—z relation therefore results in 1.5 ^ 2:1a. s ~ 2. The 
uncertainty from the t — Z relation is smaller. 

6.2. SN la Rates in Elliptical Galaxies 

For cluster ellipticals, the uncertainty in the t — Z rela- 
tion is much smaller than spirals, because [Fe/H] reaches 
~ — 1 in several 10^ years by the star burst. The parame- 
ters for SFR is constrained from the color evolution of the 
cluster ellipticals, although other star formation histories 
may be allowed for the field ellipticals. 

Like 2ia,s, the t — z relation affects zia.o- Especially a 
large uncertainty lies in Zf . In our calculations, we assume 
Zf ~ 4.5 and it results in the decrease of the SN la rate 



twice at zja 



1.6 and ~ 2.6, which respectively corre- 



spond to the hfetimes of the MS+WD and the RG+WD 
systems after Zf (Fig.^). However, the observational con- 
straint from the color evolution on Zf only requires Zf ^ 3 
for cluster ellipticals. With other cosmological models and 
Zf , the uncertainty in the t—z relation results in zia,o ~ 2 
and 1 < zia,e 2 for the MS-hWD and the RG-hWD sys- 
tems, respectively. 

For field ellipticals, there are large uncertainties in both 
t — z and t — Z relations owing to little information from 
observations, including the possibility that there is little 
difference between cluster and field ellipticals. To discuss 
quantitatively the uncertainties in zia,c in field, further ob- 
servational information would be needed. 

6.3. Total SN la Rates 

If the many of the host galaxy types are not identi- 
fied, the total rate of the spirals and elliptical compo- 
nents would depend on the parameters in more compli- 
cated manner. 

In Figure Ob, the first break of the total SN la rate ap- 



pears clearly at zia,tot ~ 1-9. This is the case if Zf ~ 4.5 
and the contribution of ellipticals (zia,o ~ 1-6) is smaller 
than spirals (zia,s ~ 1.9). If Zf ~ 3, for example, the 
break at zia,tot ~ 1.9 disappears. The second break at 
■^la.tot ~ 2.6 corresponds to the second zia,e- 

Furthermore the ambiguity in the relative mass frac- 
tions Ci of spirals and ellipticals makes a prediction of the 
first zia.tot more complicated, while the second one remains 
unchanged. If the relative contribution from ellipticals is 
larger, the first zia,tot appears at the first zia,e (e-g-, z ~ 1.6 
for Zf 4.5). Besides in our calculations the relative mass 
fractions Ci of galaxies are assumed to be constant against 
the redshift. If there is the number evolution of galaxies as 
suggested by the hierarchical clustering scenario, the star 
formation history and thus the supernova rate in ellipticals 
are similar to those for field ellipticals. 

7. CONCLUSIONS 

We have made predictions of the cosmic supernova rate 
history as a composite of the supernova rates in different 
types of galaxies. We adopt the SN la progenitor sce- 
nario including the metallicity effect which successfully 
reproduces the chemical evolution of the solar neighbor- 
hood (K98). We construct the evolution models for spi- 
ral and elliptical galaxies to meet the latest observational 
constraints such as the present gas fractions and colors 
for spirals, and the mean stellar metallicity and the color 
evolution from the present to z ^ 1 for ellipticals. 

Owing to the presence of the two kinds of the SN la 
progenitor systems (MS-I-WD and RG-fWD) with shorter 
(0.5 — 1.5 Gyr) and longer hfetimes (2 — 20 Gyr), respec- 
tively, we can explain the difference in the relative ratio of 
the SN la to SN II rate TZin/TZu between early and late 
types of spirals. Owing to the presence of the RG-I-WD 
systems with over 10 Gyr lifetime, SNe la can be seen 
even at present in ellipticals where the star formation has 
already ceased more than 10 Gyr before. 

We then construct the cosmic SFR as a composite of 
the SFRs for different types of galaxies, and predict the 
cosmic supernova rates for two different environments: 

1. In the cluster environment, the synthesized cosmic 
SFR has an excess at z <L 3 due to the early star 
burst in ellipticals and a shallower slope from the 
present to the peak at z 1.4, compared with 
Madau's plot. The predicted cosmic supernova rate 
suggests that in ellipticals SNe la can be observed 
even at high redshifts because the chemical enrich- 
ment takes place so early that the metallicity is large 
enough to produce SNe la at z 2.5. In spirals 
the SN la rate drops at z ~ 2 because of the low- 
metallicity inhibition of SNe la. 

2. In the field environment, ellipticals are assumed to 
form at such a wide range of redshifts as 1 ^ z 4. 
The synthesized cosmic SFR has a broad peak 
around z ~ 3, which is in good agreement with the 
recent sub-mm observation by Hughes et al. (1998). 
The SN la rate is expected to be significantly low 
beyond z ^ 2 because the SN la rate drops at z ~ 2 
in spirals and gradually decreases from z ^ 2 in el- 
lipticals. 



KOBAYASHI, TSUJIMOTO, AND NOMOTO 



9 



This work has been supported in part by the grant-in- 
Aid for Scientific Research (08640336) and COE research 
(07CE2002) of the Ministry of Education, Science, Cul- 
ture, and Sports in Japan. C.K. thanks to the Japan Soci- 
ety for Promotion of Science for a financial support. K.N. 
thanks to the participants in the workshop on "Type la Su- 



pernova" (Aspen Center for Physics, 13-24 June 1999) for 
informative discussion. We would like to thank I. Hachisu 
and M. Kato for providing us with their new results, B. 
Schmidt for fruitful discussions, and T. Kodama for pro- 
viding us with the database of simple stellar population 
spectra. 



APPENDIX 

B PARAMETER DEPENDENCE 
Present supernova rates 

The SN la rate prediction involves some uncertainties associated with the SN la progenitor model, which arc mostly 
reflected in the uncertainty in b parameters given in section 2.1. 6ms and 6rg are the parameters for the total fraction 
of the WDs (with the initial mass of 3 - 8Mq) which eventually produce SNe la in the MS-hWD and the RG+WD 
systems, respectively (see Eq.[|l^ in section 2.2 for the exact definition of b). In the present study, the parameter set of 
[bus = 0.02, 6rg = 0.05] is chosen to best reproduce the chemical evolution in the solar neighborhood. Other sets of b 
parameters can also be approximately consistent with the chemical evolutions, but provide different predictions of the 
present supernova rates as follows. 

a) supernova rates in spiral galaxies — - In Figure ^ we use the relative ratio of the present SN la to SN II rates TZia/Ti-u to 
avoid the Hq uncertainty included in the observed supernova rates. The adopted set of [&ms = 0.02, 6rg = 0.05] is seen to 
give also the best fit to TZia/TZu. If we adopt larger 6rg, T^ia/T^ii is larger. For [&ms = 0.04, 6rg = 0.04], TZia/TZu ~ 0.4 7 



and 0.26 for early and late type spirals, respectively, which are larger than the observed values (CappcUaro et al. 1999) 
For [6ms — 0.01, 6rg — 0.055], TZia/TZu ~ 0.28 and 0.21 for early and late type spirals, respectively, which are smaller. 
However, these sets of b parameters may be allowed in view of the possible systematic errors and the inhomogeneity of 
the samples. We need further updates of the present supernova rates to test our model. 

b) SN la rate in elliptical galaxies — - In ellipticals, the star formation has stopped more than 10 Gyr before and only 
SNe la from the RG-I-WD system occur at present (Fig.|^). Thus 6rg can be constrained from the present SN la rate in 
ellipticals. For example, if we adopt the sets of [6ms = 0.04, 6rg = 0.04] a nd [6ms — 0.055, 6rg = 0.01], the predicted SN 
la rates are by a factor of 2 higher and lower than the observed values ( [Cappellaro et al. 1999 ), respectively. However 



such differences from the observed rates are comparable to the uncertainty stemming from iJo, so that these sets cannot 
be precluded until more accurate iJo is obtained. 

Supernova rates at high and intermediate redshifts 

The supernova rates at high and intermediate redshifts can also be used to test our SN la model. In particular, whether 
the 6 parameters are different among the different types of galaxy could be seen from the evolution of the supernova rates 
with redshift. Observationally, we have found no significant differences in the binary fraction and the iron abundance 
between ellipticals and spirals. Theoretically, as noted in section 2.1, 6 may depend on the binary fraction and the iron 
abundance; especially, the metallicity dependence of b^Q is significant. Here we show the results for the cases that 6rg is 
larger only for ellipticals. 

Figure 9 shows the cosmic supernova rates for the set of [6ms — 0.05, 6rg = 0.02] in spirals and [6ms = 0.05, 6r,g = 0.10] 
in ellipticals. In this case the contribution of ellipticals to the cosmic SN la rate is as large as that of spirals. The redshifts 
where the SN la rates drop are the same as in Figures ||b and ||b. This is because 2:1a, o of ellipticals is determined by the 
lifetime of SNe la (i.e, the time-redshift relation) and not the iron abundance of the galaxy (i.e., the age-metallicity relation, 
which is effective only for spirals), while the difference in 6 affects only the latter relation. Whether the contribution of 
ellipticals is actually as large as spirals will be seen soon in the observations for the intermediate redshifts (B. Schmidt, 
private communication) 



In the above case, the present SN la rate in ellipticals is ~ 5 times larger than the observed rate (Cappellaro et al 



1999). However, the SN la rate could significantly decrease to the present, if the lowest mass of the RG companions to 
produce SNe la is not O.9M0 (with the lifetime of ~ 19 Gyr for Z = 0.02) but l.OMg (with the Hfetime of 12 Gyr for 
Z = 0.02), because the star formation in ellipticals has stopped 14 Gyr before in our model. Whether the longest lifetime 
of SNe la is actually younger than the age of ellipticals can be seen from the intermediate redshift observations. 
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Fig. 1. — The distribution functions of the companion mass (upper panel) and the companion lifetime (lower panel). In the lower panel, 
the solid and dashed lines are for Z = 0.002 and Z = 0.02 (solar), respectively. MS+WD and RG+WD denote that the companions of the 
white dwarfs (WD) is the somewhat evolved near main-sequence (MS) star and the red-giant (RG), respectively. 
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Fig. 2. — Star formation rates (SFR: top panel), gas fractions (middle panel), and B — V colors (bottom panel) as a function of time for 
four types of spirals : SOa-Sa (solid line), Sab-Sb (long-dashes line), Sbc-Sc (short-dashed line), and Scd-Sd (dotted line). The crosses show 
the present observed values of colors and gas fractions. The present B-V colors are taken from Roberts & Haynes (1994). We use the gas 
(i.e., HI+H2) fractions which are normalized by the present blue luminosity of the galaxy to avoid the uncertainty in the fractinna of the dark 
matter. The HI mass is taken from Roberts & Haynes (1994), and the H2 mass is derived from the H2/HI ratios ([Oasoli et al. 1998). 
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Fig. 3. — (a) TViq pggc^^ro fn}r,r <nyr.1iitir.Ti prorlir-torl h^r the niodel of cluster elUpticals from the present to 2 ~ 1, compared with the 
observational data ( Stanford, Eiscnhardt fc Dickinson 199q ). (b) The distribution function of the formation r nnrh adopted in thp mn dpl of 
field ellipticals. The observational data arc estimated from the spectra of ellipticals in the Hubble Deep Field ( Franceschini et al. 199^). 
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Fig. 4. — The ratio of the SN la rate to SN II rate TZiafR-ii as a function of time for spiral galaxies. The upper and lower lines show 
the results for the early-type spirals SOa-Sb and for late-type spirals Shr-f^H rpgpertiwly The solid, dashed, arid HnttpH linea are ralmlaterl 
with our SN la model, the single delay-time model with tia ~ 1.5 Gyr (Yoshii et al. 1996), and the DD model (Tutukov & Yungelson 1994), 
respectively. The observational data are taken from Cappellaro et al. (1999). 




Fig. 5. — The SN la rate as a function of time in elliptical galaxies. Lines are the same as in Figure ^. The observational data is taken 
from Cappellaro et al. (1999). 
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square), (c) The cosmic supernova rate per volume (solid line) as a composite of those in 
spirals (long-dashed line) and ellipticals (short-dashed line). The observational data for the cosmic SN la rate is taken from Pain (1999). 
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Fig. 7. — (a) The mgmir SFR alnng; rpH sViift which is a connected line of the observed data. The sohd Hne takes into account the correction 
of the dust extinction (Pettini et al. 199^) and the dotted hne does not. (b) The iron abundance in the gas of spirals (long-dashed line) and 
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Fig. 8. — The same as Figure nl but for ellipticals whose formation epochs span over 1 ^ 2 ^ 4; this might correspond to field ellipticals. 




Fig. 9. — The cosmic supernova rate along redshift in clusters^a) and field (b) with [feMS = 0.05, fenG = 0.02] for spirals and [bjvis = 
0.05, ftjiG = 0.10] for ellipticals. Symbols are the same as in Figure tfo. 
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Table 1 

The input parameters in equations (^) and (^). 





Tb [Gyr] 


n [Gyr] 


(1) 


(2) 


(3) 


E 


0.1 


0.1 


SOa, Sa 


2.32 


2.35 


Sab, Sb 


2.19 


5.61 


Sbc, Sc 


2.39 


10.92 


Scd, Sd 


2.19 


60.5 



Col.(l). — Galaxy type. 

Col. (2) (3). — Timescales of the 
star formation and inflow in Gyr. 



Table 2 

The calculated quantities for spirals and ellipticals at 15 Gyr. 





/g 




[M/H]g 


[M/H]s 


[Fe/H]g 


[Fo/H]b 


[Mg/Fe]g 


[Mg/Fe]B 




nu 




TZu 




B-V 


U ~B 


M/Lb 


M/Lv 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 




(10) 


(11) 


(12) 


(13) 


(14) 


(15) 


(16) 


(17) 


E 


0.095 


0.905 


0.176 


-0.212 


0.831 


-0.411 


-0.777 


0.339 




0.000 


0.007 


0.000 


0.084 


0.924 


0.339 


11.323 


8.799 


SOa, Sa 


0.031 


0.968 


0.323 


-0.160 


0.607 


-0.026 


-0.215 


-0.032 




0.072 


0.034 


0.425 


0.198 


0.776 


0.106 


5.880 


5.238 


Sab, Sb 


0.070 


0.862 


0.082 


-0.184 


0.315 


-0.050 


-0.151 


-0.032 




0.174 


0.051 


0.712 


0.209 


0.631 


-0.056 


4.105 


4.178 


Sbc, Sc 


0.098 


0.648 


-0.033 


-0.218 


0.164 


-0.092 


-0.109 


-0.023 




0.225 


0.053 


0.853 


0.202 


0.541 


-0.137 


3.796 


4.199 


Scd, Sd 


0.038 


0.182 


-0.124 


-0.235 


0.051 


-0.112 


-0.081 


-0.019 




0.095 


0.019 


0.953 


0.194 


0.472 


-0.189 


10.085 


11.880 



Col.(l). — Galaxy types. 

Col.(2)(3). — Gas and steUar fractions. For eUipticals, /g includes the gas ejected in the galactic 
wind. 

Col.(4)(6)(8). — Gas metallicities [M/H]g = logZ/ZQ, iron abundances of gas [Fe/H]g, and mag- 
nesium to iron ratios of gas [Mg/Fe]g. 

Col.(5)(7)(9). — Mean stellar metallicities [M/H]s = log Zb/Zq, mean stellar iron abundances 
[Fe/H]a, and mean stellar magnesium to iron ratios [Mg/Fe]s. 

Col.(10)(ll). — SN II and la rates per mass in the unit of [/Gyr/lO^Afo]. 

Col.(12)(13). — SN II and la rates per luminosity in the unit of [/Century/IO^^LBQ]. 

Col.(14)(15).— B a.ndU - B colors. 

Col.(16)(17). — Mass-to light ratios in B- and V-band. 



